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Radiative decays of X(3872) with J PC = 1 ++ are studied in the coupled-channel 
approach, where the cc states are described by relativistic string Hamiltonian, while 
9-; for the decay channels DD* a string breaking mechanism is used. Within this method 

a sharp peak and correct mass shift of the 2 3 Pi charmonium state just to the D°D*° 



(N 



(N 



X 



threshold was already obtained for a prescribed channel coupling to the DD* decay 
channels. For the same value of coupling the normalized wave function (w.f.) of 



> 

qq , X(3872) acquires admixture of the 1 Pi component with the weight c\ = 0.153 

■>sj- ' 8.8°), which increases the transition rate T(X(3872) — > J/ip'y) up to 50-70 keV, 

making the ratio R = Q^x{^72)^t/^\) = ^-8 ^ 0-2 (th) significantly smaller, as 



compared to R ~ 5 for X(3872) as a purely 2 3 Pi state. 



I. INTRODUCTION 

The AT (3872) was discovered by Belle as a narrow peak in J/ipirTt invariant mass dis- 
tribution in decays B — > J/i/jttttK [1] and later confirmed by the CDF, DO, and BaBar 
Collaborations [2]. It has several exotic properties, very small width T < 2.3 MeV and the 
mass very close to the D°D*° threshold [3], [4]. The even charge parity C = + of A(3872) 
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is now well established [5], while two most plausible assignments for its quantum numbers, 
jpc = 1+ + and 2 -+ ) are st m discussed [6], [7]. 

To understand the nature of X(3872) a special role belongs to radiative decays, 
X(3872) and X(3872) ^'7. The first evidence for the decay X(3872) -> 

was obtained by Belle [8] and confirmed by BaBar [9]; later the BaBar has observed radiative 
decay X(3872) -> ^'7 with the branching fraction ratio P = = 3A± 1A t 10 ]- 

Knowledge of this ratio is of special importance for theory, because the rates of these ra- 
diative decays vary widely in different theoretical models [11]-[15]. If X(3872) is considered 
as a conventional 2 3 P 1 charmonium state, then the characteristic value of R is rather large, 
R ~ 4 — 6 [12]- [14], being in general in agreement with the BaBar number 3.4 ± 1.4. In 
molecular picture the radiative decay X(3872) — > ip'j is suppressed and the ratio R should 
be much smaller [11]. 

However, in 2010 on a larger sample of decays B — > X(3872)K) the Belle has not found 
evidence for the radiative decay X(3872) — > ^'7, giving the upper limit R < 2.1 [16]. This 
number does not agree with the representation of X(3872) as a purely 2 3 Pi charmonium 
state. 

Existing experimental uncertainty calls for new studies of coupled-channel (CC) effects 
for X(3872). In [13] the authors have anticipated "... a significant DD* component in 
X(3872), even if it is dominantly a cc state". For that the 3 P model was used in [13] and 
the Cornell many-channel model was considered in [14]. However, in spite of these many- 
channel calculations there predicted values of R have appeared to be close to those obtained 
for X(3872) as a purely 2 3 Pi charmonium state: R = 5.8 in [13] and R = 5.0 in [14]. 

Here we consider X(3872) with J PC = 1 ++ in the CC approach, where a coupling to the 
DD* channels is defined by the relativistic string-breaking mechanism, which was already 
applied to X(3872) in [17], [18], explaining it as 2 3 Px charmonium state shifted down and 
appearing as a sharp peak just at the D°D*° threshold. Besides, the scattering amplitude 
and a production cross section were calculated there, being in qualitative agreement with 
experiment. Here we apply this method for calculations of the radiative decay rates for 
X(3872) and show that due to the same CC mechanism (with the coupling of the same 
strength) an admixture of the l 3 Px component to the X(3872) w.f. appears to be not 
large, ~ 15%; nevertheless, this component strongly affects the value of the partial width 
T 1 = T(X(3872) -> J/ijrf) and decreases the ratio R. 
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II. COUPLED-CHANNEL MECHANISM 

We use here the string decay Lagrangian of the 3 P type for the decay cc — > (cq) (cq) [18] : 



■'sd 



$ q M u i/> q d 4 x (1) 



where the light quark bispinors are treated in the limit of large m c , which allows us to go 
over to the reduced (2 x 2) form of the decay matrix elements (m.e.). Also to simplify 
calculations the actual w.f. of cc states, calculated in [19] with the use of the relativistic 
string Hamiltonian (RSH), is fitted here by five (or three) oscillator w.f. (SHO), while the D 
meson w.f. is described by a single SHO term with few percent accuracy with the parameter 
j3 ~ 0.48. In this case the factor in (1) is M u ~ ^ ~ 0.8 GeV, which produces correct 
total width of -?/>(3770) and it will be used below. 

The transition m.e. for the decays (cc) n — > (DD), (DD*), (D*D*) are denoted here as 
n — > rii-, n 3) an d in the 2x2 formalism this m.e. reduces to 

7 f <i 3 q 

^nn 2 n 3 (p) = ^= J J/l23^3*n(P + (q)^n3 (o) ■ ( 2 ) 

Here 7 = <m + ™y D+eo> , where average of the Dirac denominator (with scalar confining 
potential U = ar and vector potential Vd = — f^) is calculated and yields 7 = 1.4. The 
factor yi23 contains a trace of spin-angular variables (for details see [18], [20]). 

The intermediate decay channel, like DD*, induces an additional interaction "potential" 
Vcc{q,q',E) (here the quotation marks imply nonlocality and energy dependence of this 
potential): 



where 



(3) 



^n 2 n 3 (q,p) = -^=J/I23(q,p)^n a (q-P)^n 3 (q-P)- ( 4 ) 

Using (3) and (4) one can find how the energy eigenvalues (e.v.) and the w.f. of a state 
(cc) n ) change due to the interaction Vcc- in particular, in the first order of perturbation 
theory one has 
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=E n + w nn (E n ), (5) 

> = ^» + £ 't^T^- (6) 
where £■„ refer to the unperturbed (cc) n system and the m.e. w nm (E) is 

<i 3 q d 3 q' 



^) = / ^(q)^c(q,q'^)^(q')^^ = (7) 



d^P Jnnmz (f*) ^mnmz (p) 



( 27r ) 3 6r 3 E-E n2m:i ( P ) 

Notice, that if the CC interaction is strong, then one should take into account this 
interaction to all orders, summing the infinite series over Vcc (or w nm ). As a result, one 
obtains the full Green's function for an arbitrary QQ system (our result formally coincide 
with those from [21], although differences occur in concrete expressions for w nm , because 
our interaction and decay mechanism differ from those in [21]): 

G QQ (l, 2;E) = Y: <Ol)(£ -E + t&U*$°(2), (8) 

n,m 

here E nm = E n 5 nm . Then the energy e.v. are to be found from the zeros of the determinant 

det(E - E -w) =0. (9) 

If the mixing of the state n with other states, m ^ n, is neglected, one obtains a nonlinear 
equation for the e.v. E* with a correction: 

E^=E n + w nn (E^). (10) 

Note that E* can be a complex number and occur on the second sheet in the complex 
plane with a cut from the threshold n^n-i to infinity. The expression (10) was used in [17] to 
find the position of the pole for shifted 2 3 P\ charmonium state, while its unperturbed value 
E 2 was calculated with the use of RSH [22] to be E 2 = 3948 ± 10 MeV. In [17] the decays 
have included D Dq and D + D~ channels and 7 was used as a free parameter. Then the 
position of the resulting pole (10) and the production cross section have been calculated, 
giving the pole position exactly at the D Dq threshold for 7 = 1.2, which is close to expected 
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number 7 = ^ + ^_y D+eo ^ = 1.4 for M w = 0.8 GeV. (From Fig. 2 in [17] one can see that for 
this value of 7 the production curve agrees well qualitatively with experiment.) 

Knowing the mass shift of X(3872), one has also to find next order corrections to the 
w.f., thus going beyond no-mixing approximation. At first, we use perturbation theory, 
when admixture of the states m 3 P 1 = m\ ++ to the 2 3 P 1 state is given by the (6), 

cS> = -=^, m = 1,3,4,.... (11) 

Calculations give w 2 i = 0.085 GeV, w 23 = -0.0008 GeV, while E 2 - E 1 = 0.425 GeV, so 
that main contribution comes from the m — 1 (1 3 Pi) state with the mixing parameter 

c? = " °- 2 °- W 

This correction S^ 1 is not small and calls for a more accurate calculations, beyond pertur- 
bation theory. 

To this end we first write general expressions for the yield of particles 7, 7r, p, u etc. from 
the system, originally born as a (QQ) system in e + e~ or B meson decay by an operator B. 
The Green's function for the system can be written as 

offi = L ' ) (**-)• < 13 > 

where e.g. (Bip n ) ~ VVi(0) for the e + e~ production. We now can attribute particle (i) pro- 
duction from Q, Q lines adding the corresponding self energy parts to E, (E) nm = E n 5 nm + 
Ylnm(E), and the production from light quark lines to w mn (e), w mn (E) — > w mn (E)+Wm n (E). 
Then the yield for particles % can be written as 

Y t (E) = £ m n ) ( 1 ) ^(E« + ( * — ) ^ (fl^)*- (14) 
\E-E + wJ nm 2t \ E -E + wJ lq 

One can define now effective (QQ) wave function ^qq, which actually participates at the 
vertex of emission of particles i, 

fc,i \& — & + w / kl k 
Keeping only two eigenfunctions, one has 

ax = ^ Wr) - ^(^2) (16) 
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E x - E + w u - w 12 ~ 

a 2 = ^ (Brh) ~ ^ WO, 



(17) 



where det = det(E — E + w). 

^From (21) one obtains the ratio 



W 2 l + (E R - E) 



(%) 
(BV2) 



ai 



° 2 ^"-n + ^gg 



Ai = - E x = 362 MeV. 



(18) 



Neglecting and for E = E R = 3872 MeV one obtains 



w 2 i 



a 2 Ai - tun 



0.179, 



(19) 



while approximating in £?-decay production 



(JNn) „ ^1(0) 
(BV2) ~ V>a(0) 



0.85 one obtains 



— = 0.155. 

a 2 



(20) 



One can see in Table I, that the n = 3 3 Pi state gives a negligible admixture. The values 
of w mn for E = En are computed according to Eq. 7 with w.f. obtained in [19] and are 
given in Table I. 

TABLE I: The m.e. w nm (in GeV) between n 3 P\ and m 3 P\ states for two approximations of exact 
w.f. 



nm 


11 


12 


22 


32 


("'nm 


-0.320 


0.122 


- 0.099 


-0.0003 


5SHO 












-0.319 


0.121 


- 0.098 


-0.0011 


3SHO 











Then using (20) the w.f. of X(3872) can be presented with a good accuracy as 

y?(X(3872)) = 0.988 ip(2 3 P ± ) + 0.153 </?(! 3 Pi) (21) 



III. RADIATIVE DECAYS 



Electric dipole transitions between an initial state (i) n 3 Pi state and a final (f) state 
m 3 Si are defined by the partial width [7], [24], 
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T( l ^ 1 + f) = ^ae 2 Q E 3 (2J f + l)Sf f \S tf \ 2 , (22) 



where the statistical factor = is 




Sf f = max (/,/')< > • (23) 



For the transitions between the n 3 Pj and m 3 Si (m 3 Di) states with the same spin S = 1, 
the coefficient Sfj = |(^)- 

To calculate m.e. Sif we use RSH H [22], which is simplified in case of heavy quarkonia 
when one can neglect a string and self-energy corrections, arriving at a simple form (widely 
used in relativistic potential models with the constituent quark masses in the kinetic term 
[25],[26]): 

H = 2 v /p 2 + ml + V B (r). (24) 

By derivation, in (24) the mass of the c quark cannot be chosen arbitrarily and must 
be equal to the pole mass of a c quark, m c ~ 1.42 GeV. The pole mass takes into ac- 
count perturbative in a s {m c ) corrections and corresponds to the conventional current mass 
m c (m c ) = 1.22 GeV [27] (here m c = 1424 MeV is used). 

The potential V B (r) taken, 

V B (r) = or - *|M (25) 

contains the string tension (er = 0.18 GeV 2 ), which cannot be considered as a fitting param- 
eter, because it is fixed by the slope of the Regge trajectories for light mesons. In the vector 
strong coupling the asymptotic freedom behavior is taken into account with the QCD 

constant A B , which is defined by A^: A B (n f = 4) = 1.4238 A S 5(n / = 4) = 370 MeV, the 
latter is supposed to be known; in our choice A B (nf = 4) corresponds to Kj^(nf = 4) = 
261 MeV. At large distances a B (r) freezes at the value a cr it = 0.60. 

Then for a given multiplet nl the centroid mass M cog (nl) is equal to the e.v. of the 
spinless Salpeter equation (SSE): 



H ipni = M (nl)ip nl . 



(26) 
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We have calculated M cog (nl) in two cases: in single-channel approximation, when 
M cog (2P) = 3954 MeV was obtained, and also taking into account creation of virtual loops 
qq, which are important for the states above the open charm threshold and give rise to flatten- 
ing of confining potential [28]; in the last case M cog (3P) = 4295 MeV, M cog (2P) = 3943 MeV 
(which is by 9 MeV smaller than without flattening effect), and then due to fine structure 
(FS) splittings the mass M(2 3 P 1 ) = 3934 MeV is calculated. 

For a multiplet nP a spin-orbit a so (nP) and tensor t(nP) splittings are calculated here 
taking spin-orbit and tensor potentials as for one-gluon-exchange interaction, although as 
shown in [29], second order (aj s (fj,)) corrections appear to be not small for the IP multiplet; 
their contribution can reach < 30% . 

a so (nP) = ^a fs (r~ 3 ) nP - + t(nP), (27) 

t(nP) = ^{r- 3 ) nP . (28) 
o uj c 

We take here atf s — 0.37, which provides precise description of the fine-structure (FS) 
splittings for the 1 3 Pi charmonium multiplet, if second order corrections are taken into 
account [29]. For the 2P multiplet the masses: M(2 3 P 2 ) = 3963 MeV, M(2 3 Pl) = 3934 
MeV, M(2 3 P ) = 3885 MeV, M(2 1 P 1 ) = 3943 MeV are obtained. Notice, that one cannot 
exclude that for the states above open charm threshold the FS splittings may be smaller or 
totally screened due to coupling to the DD* channel, and even the order of the states with 
different J may be changed. 

Below we use the following mass differences: 

M cog (2P) - M cog (lP) = 425 MeV, M co ,(3P) - M cog (lP) = 770 MeV, 

M cog (3P) - M cog (2P) = 350 MeV, (29) 

In Table 2 the m.e. Sif (in GeV) between n 3 Pi {n = 1,2) and m 3 Si states are given; in 
some cases, if the value of Sif is small and results strongly depend on a/ s (/i) used, we give 
two variants: first, with "normal" FS splittings and in second case FS effects are totally 
suppressed. 

To control an accuracy of our calculations in Table 2 we give also the partial widths of 
the dipole transitions: 1 3 Pi ->■ 2 3 S'i ->■ Xd(3510), and 1 3 D t ->■ x cl (3510), and in 
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TABLE II: El transition rates. The m.e. (X(3872)|r|n 3 Si) (n = 1,2) includes admixture from 
the l 3 Pi component with c\ = 0.153; experimental data from [30]-[32]. 



Transition 




Sif 


Sif 


r( 


i ->• /) (keV) 


. El 


/ 


(MeV) 




(GeV)- 1 


this paper 


BG[13] 


experiment 


l 3 Pi(3510) 


l 3 Si(3097) 


389 


l 

9 


1.927 


315 




317±25[30] 


2 3 Pi(3872) 


l 3 5i(3097) 


697 


1 

9 


0.104°) 
0.216 6 ) 


5.3 
22.8 


11.0 




V(3872) 


l 3 5i(3097) 


697 


1 

9 


0.396 c ) 


76.6 












> 0.292 d ) 


> 41.7 






2 3 S'i(3686) 


l 3 Pi(3510) 


171 


1 

9 


-2.104 


31.9 




30.6 ± 2.2[31] 


2 3 Pi(3872) 


2 3 Si(3686) 


181.5 


1 

9 


3.02 


78.6 


63.9 




V(3872) 


2 3 S 1 (3686) 


181.5 


1 

9 


2.70 


62.8 






l 3 L>i(3770) 


l 3 Pi(3510) 


252.9 


1 

18 


2.767 


89 


199 


70 ± 17 [32] 
80 ± 24 [27] 


2 3 P 2 (3872) 


l 3 Z?i(3770) 


97.7 


1 

18 


-2.776 


5.2 


3.7 




V(3872) 


1 3 Z?!(3770) 


97.7 


1 

18 


-2.32 


3.6 







a ) otf s = 0.37 in spin-orbit potential. 
fc ) FS interaction is totally suppressed. 

c ) Both admixture of the l 3 Pi state and FS splittings with af s = 0.37 are taken into account. 

d ) The lower limit refers to the case when FS interaction is totally suppressed. 

all three decays good agreement with existing experimental data [30]- [32] is obtained (here 
the 2 3 Si and 1 3 Di states are identified with ip' and ip"). In our calculations the 2S — ID 
mixing in not taken into account. 

lYiom. Table 2 one can see that for X(3872) as a purely 2 3 Pi state, the transition rate 
T 1 = r(2 3 P x — > J/V't) strongly depends on FS potential used. For this transition the 
m.e. S 2 i = 0.10 GeV" 1 is small, if in spin-orbit potential "normal" otf s — 0.37 is used, and 
Ti = 5.3 keV is also small, being less than in [13], where a larger £ 2 i = 0.15 GeV" 1 and 
Ti = 11 keV were calculated. 

In the case when FS interaction is neglected, or suppressed, then m.e. £ 2 i = 0.216 GeV -1 
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is larger (in [13] £ 2 i ~ 0.276 GeV -1 ), and the partial width T 1 = 22.8 keV is 4 times larger. 
In both cases discussed a contribution from the 1 3 Pi component was also neglected and the 
ratio R is large, R > 4. Notice, that in [33], using an analogy with the radiative decays of 
Xbj{2P) in bottomonium, a smaller value of this ratio, R = 1.64 ± 0.25, was predicted. 

However, if admixture from the 1 3 Pi state, as in (21), is taken into account, then the 
transition rate F 1 increases, independently of a strength of FS interaction used. With 
Ci = 0.153 and suppressed FS interaction we obtain the lower limit for r 1; 



This transition rate reaches a larger value, T = 76.6 keV, if in FS potential the same aj s = 
0.37, as for the IP states, is used. On the contrary, the partial width T 2 = T(X(3872) — y 
ip ,r y) decreases (by 20%), owing to admixture ci in the w.f. of X(3872) and negative m.e. 
(2 3 S'i|r|l 3 Pi). As a whole, the ratio R is becoming smaller and can change in wide range: 



where the upper limit refers to the case when FS interaction is strong and theoretical error 
comes from a variation of aj s , while the lower limit refers to the case when FS interaction 
is suppressed. 

These values of R are in agreement with the Belle measurements of the X(3872) radiative 
decays where a restriction R < 2.1 was observed [16]. At the same time our limit, R < 
(0.8 ± 0.2), does not agree with R = 3.4 ± 1.4 from the BaBar experiment, being also much 
smaller than in many-channel calculations [13], [14], where the ratio R ~ 5 was obtained. 

Thus we conclude that precise measurements of R are of great importance for under- 
standing the nature of X(3872): firstly, for definition of admixture of the l 3 Pi component 
in the w.f. of X(3872) and secondly, for understanding of FS effects in higher resonances, 
which lie above open-charm threshold. 

Notice, that the partial width Ti(2 3 Pi(3872) — > J/ip^) is not very small even for a pure 
2 3 Pi state, if a contribution from spin-orbit interaction is suppressed, and the partial width 
increases 4 times: 



r!(X(3872) ->■ J/^7) > 41.7 keV. 



(30) 



0.53 < R< (0.8 ±0.2) (th), 



(31) 



ri 0) (2 3 P 1 (3872) -> Jji)i) = 22.8 keV, 



(32) 
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while in this case rj> ) (2 3 P 1 (3872) — > ip'j) = 85.5 keV increases only by ~ 8%; however, 
their ratio remains large, 

R (2 3 P 1 ) = 3.75. (33) 

The situation changes, if FS interaction is suppressed, but the w.f. of X(3872) contains 
admixture of the 1 3 P 1 state (21). Then the m.e. (J*T(3872)|r| = 0.512 GeV-1 and 

Ti = 128 keV reaches the maximum value, but the m.e. (X(3872)|r|-?//7) = 2.80 GeV 
changes by only ~ 10% and T 2 = 67.6 keV; so that their ratio has a minimal value: 

Rmin) = 0.53. (34) 
To define FS splittings of the nP multiplets we use following m..e. 

(2P|r _1 |lP) = 0.134 GeV, (2P\ r - 3 \lP) = 0.123 GeV 3 , 

(3P]r _1 |lP) = 0.080 GeV, (3P\r~ 3 \lP) = 0.112 GeV 3 , (35) 

(3P|r _1 |2P> = 0.133 GeV, (3P\r- 3 \2P) = 0.134 GeV 3 . 
Then for ctf s — 0.37 one finds 

(X(3872)|r|J/^ 7 ) = 0.104 GeV^ 1 , (36) 

which is two times smaller than in the spin-average case, when (X(3872)|r| J/ip'y) = 0.216 
GeV -1 , and the partial width is large, 

rx(X(3872) -> J/rjn) = 76.6 keV. (37) 
The m.e. (X (3872)\r\ip'^) = 2.70 GeV" 1 in this case, giving 

r 2 (X(3872) ->■ V't) = 62.8 keV, (38) 

and the ratio R = 0.82. However, if a stronger af s > 0.45 is used, this ratio may reach a 
larger value, ~ 1.0, so that R can vary in the range, 



R = 0.8 ±0.2 (th). 



(39) 
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Thus, while both spin-orbit splitting and admixture C\ from the 1 3 Pi state are present, then 
the partial widths T 2 and F 1 turn out to be of the same order. 

In our calculations above we have disregarded the contribution of the 7 emission from 
the DD* intermediate state in the radiative decays of X(3872) into J/ip or if)'. To estimate 
this part of the 7 emission we refer to calculations done in [34]. It was found there that 
the channel DD* contributes to the J/ip'y final state less than 3.6 keV and to the ^'7 final 
state less than 0.01 keV, i.e. these contributions are smaller as compared to changes in 
corresponding partial widths due to variations of the coupling aj s in the range 0.25 — 0.45 
(see Table 2). Therefore in this paper we have disregarded possible effect of 7 emission from 
DD* intermediate states. 

IV. CONCLUSIONS 

We study the exotic charmonium state X(3872) with J PC = 1 ++ in the CC approach, 
where a coupling to the DD* channels is determined by the parameter-free string-breaking 
mechanism. Due to this coupling the 2 3 Pi charmonium state is shifted down to the D°D* 
threshold and its w.f. acquires admixture from the 1 3 Pi cc state. Such mixing of the 2 3 Pi 
and l 3 Pi states is not large, corresponding to the mixing angle = 8.8°. 

Owing to this admixture the transition rate ri(X(3872) — > J/ip'y) increases several times 
and reaches the value in the range 45 — 80 keV. At the same time the transition rate 
r 2 (X(3872) — > ip'ry) decreases by ~ 15%. As a result their ratio has following features: 

1. The ratio R = 0.53, if spin-orbit interaction is totally suppressed. 

2. The ratio R = 0.82, if spin-orbit interaction is defined by the FS coupling, aj s ~ 0.37 
and can reach the larger value ~ 1.1 for a larger af s . 

The partial width of X(3872) — > ^"7 appears to be small, ~ 4 keV. 

Our calculations support the Belle result that R(exp.) < 2.1, while a larger number, 
R = 3.4 ± 1.4 puts an additional restrictions on the value of admixture c\ in the w.f. of 
X(3872). 
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